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Q-i! Abstract 



Simple approximate expressions for the relative flavor fluxes of energetic cosmic neutri- 
nos detected on Earth, are presented in terms of their initial fluxes at the surface of the 
producing cosmic sites, assuming the neutrino mixing angles to lie in the experimentally 
favored region. These expressions highlight the main sensitivities to the initial produc- 
tion fluxes as well as to small variations of the mixing angles within the experimentally 
preferred region, thus providing simple methods to disentangle these physical quantities 
from cosmic neutrino data. This is more so, due to some striking mathematical properties 
of the relative neutrino flavors at (S23 = 1/V2, S13 = 0), which somehow characterize the 
whole experimentally preferred region. We also assess the quality of our approximate 
expressions through a numerical comparison with the exact results. 
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The relative flavors of energetic Neutrinos reaching the Earth, after they have been emitted 
at various cosmic sites, provide useful information on the physical conditions there. Such 
very energetic neutrinos, approaching the 10 3 TeV [T], or even the 10 6 TeV scale [2JE1II1E], 
may be generated in extra-galactic sites like Gamma Ray Bursts and Active Galactic 
Nuclei (AGN). Galactic candidates that may emit neutrinos of up to 100 TeV have also 
been identified at distances of at least 2.6kpc [HIE]- [Exploding galactic Supernovae may 
also induce observable neutrino fluxes with energies in the few MeV range [7j.] 

It is commonly believed that these neutrinos are produced mainly through the decay 
of high energy and K (and may be also some D) mesons, which implies that the initial 
relative neutrino flavors at the cosmic sites satisfy 

F e ° = 1/3 , F® = 2/3 , F° = , (1) 

[HI HI El 121 IE] • This case is referred to below, as the canonical case. It may be useful 
to remember though, that our present understanding of the mechanisms for generating 
high energy neutrinos is rather primitive, and sites may exist in the Universe where the 
produced neutrinos have a different initial structure [21 E] ■ Therefore, the measurement of 
the relative intensities of the various neutrino flavors on Earth, may provide useful direct 
information on the mechanism responsible for their generation in the Cosmos and may 
possibly lead to the discovery of some kind of New Physics [TP] . 

Once the various neutrino flavors produced inside some cosmic object reach its surface, 
they propagate oscillating through space, as dictated by the vacuum oscillation formal- 
ism 1 . Some of the intrinsic properties of the flavor oscillations are, however, not very easy 
to readout from the general vacuum oscillation formulae, especially when relying only on 
numerical scans. These properties can be of importance in the perspective of disentan- 
gling, through the measurement of relative fluxes on earth, the astrophysical uncertainties 
encoded in the initial flavor fluxes on the cosmic sites, from the particle physics features 
in the neutrino sector. It will therefore be useful to have simple formulae giving the ob- 
servable relative neutrino flavors F e , F^, F T on Earth, in terms of the initially produced 
ones F e °, F° F° at the surface of the cosmic object. 

The aim of the present paper is to give such simple analytical expressions, assuming 
only three active neutrino flavors which propagate oscillating among themselves [TTj : 
(in particular no cosmic neutrino decay is assumed J2])- These expressions have the 
advantage of identifying specific sensitivities and degeneracies in the mixing angles, and of 
encoding some essential qualitative and quantitative properties of the full formalism. Some 
of these properties turn out to be highly non-generic consequences of the concomitance 
of the physically favored values for the mixing angles s 23 , s 13 , and the initial fluxes. 

To derive the aforementioned formulae, we take into account the basic experimental 
characteristics of the neutrino masses and mixings. These are summarized as follows: 
The recent SNO ^3] data combined with those of Super-Kamiokande Jl] strongly favor 



the LMA MSW [HI] solar solution with three active neutrinos and 8 12 ~ tt/5.1 and 
1 For a review see e.g. 



2 



2 2 

mj — m 



2 | ~ 5 x 10 5 el^ 2 [T^j- The atmospheric neutrino [17j data imply # 23 — 7r /4 and 
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2.5 x 10~ 3 eV 2 ; while the CHOOZ experiment constrains 13 < 0.1, H3 El- 
Defining Sjj = sin (9jj and choosing the "central values" 

1 1 

one has 

S12 = Sj 2 + 5si2 , S 23 = S 23 + <5s 2 3 J S 13 = «13 + ^13 (3) 

where [11111111111112], 

-O.H <^ 12 < 0.04 , - 0.12 <5s 23 < 0.10 , (4) 
0<^ 13 <0.1 , (5) 

For realistic neutrino mass differences, and neutrino energies in the range E < 10 6 TeV, 
the vacuum oscillation lengths Ay = 47rE/\mj — m||, always satisfy Ay < lpc, which is 
much smaller than the distances to all cosmic neutrino emitting sites, beyond our solar 
system (Hj. Consequently, the number of oscillations performed by the cosmic neutrinos 
before arriving at the Earth, is so large, that sin 2 (7rL/Ajj) averages to 1/2, and the CP- 
violating contributions vanish. 

Studying the properties of the relative neutrino fluxes in the s 23 ,si3-piane, for any 
fixed values of si 2 and cos 5, we remark that the point (s 23 , Si 3 ) = (l/v^2, 0), (lying within 
the experimental range of Eqs.(j21 El EJ), presents the striking property of being the 
unique point of this s 2 3Si 3 -plane 2 , where a common (s 12 - and cos5-dependent) direction 
exists, along which all three fluxes F e , F^, F T are stationary. 

As an example, we expand the standard vacuum oscillation formulae to first order in 
S13 and (Ssi2, Ss 23 ) defined in ©, getting the relative neutrino fluxes on Earth, 
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(6) 



\/2 

Ss 123 = — {k5s 23 - 5s 13 cos5) , (7) 
y 



with k = 4, and (_F e , F®, F®) being the initial neutrino relative flavors at the cosmic site. 
The experimental constraints (jH EJ) then imply 

- 0.09 < 5s 123 < +0.07 . (8) 



2 At least for s 2 3 < 0.5. 
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In writing Eq.fjHJ) we took into account the unitarity relation 

F e + F. + F T = F° e + F° + F T ° = 1 



(9) 



where the right hand side is just a normalization. 

The set of equations (0) clearly indicates that for s\ 2 = l/v 7 ^, all three relative fluxes 
on Earth are stationary along the direction k[s 23 — 1/V2] — s 13 cos 5 = passing through 
the point ($23, S13) = (l/v2, 0) of the S23<si3-plane 3 . An alternative way of expressing this 
"stationary along a direction" property is to observe that 5s 2 3 and S13 enter Eqs.dHJ) only 
through one specific combination, like 5si23 of (J7J). This signals an approximate degeneracy 
in the sensitivity to these mixing angles, which would affect their reconstruction from 
experimental data. 

If second order effects are retained, then dependencies on all four mixing angle com- 
binations 5si2, 5s23, <5si23 and S13 appear, leading to 

,(3F e °-l) 

'13 J" 
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I26sj 23 + 5s 2 3(—5s 2 3 - 2\fQ8si 2 ) 
y 



Ss 123 (-V35s 12 - 12V25s 2 3 - 1)J (3F e ° - 1) 



8(3<^ 23 - 3^5123^23 + 254) (3F° - l ) 



(10) 



We next turn to the discussion of three interesting specific cases. 

Equipartition condition for fluxes on Earth. There is in fact still another reason that 

makes the specific point (s 2 3, S13) = (l/v2, 0) unique in the S23Si3-plane. To explain this, 
we study in general terms the conditions required for the initial fluxes and the mixing 
angles, in order that flavor fluxes on Earth become equal; i.e. F e = F^ = F T , a property 
referred to as flux equipartition. It is indeed noteworthy that, [besides the Supernovae 
case where flux equipartition on Earth is an immediate consequence of the (approximate) 
initial flux equipartition 4 ] , the canonical case also leads to equal neutrino fluxes on Earth, 
provided the neutrino mixing angles are appropriate. 

3 For a given central value sf 2 , the stationary direction is determined by k = (2\[2s\- 1 \/\ — (sj 2 ) 2 ) /(1 — 
2(sJ 2 ) 2 ), which indeed leads to n = 4 for s\ 2 = l/y/3. 
4 This case is discussed in more detail below. 
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Using the exact results for large distance oscillations, one finds that 



2 l-F°-2i? 2/3-(F° + F°) 

F o_ F o + (F r 0-F0)(l- s y ' [ } 

is a necessary condition 5 for equipartition of the neutrino fluxes on Earth. This condi- 
tion shows clearly that equipartition requires a specific correlation between two sets of 
physically independent quantities, namely the mixing angles S23, S13 and the initial flavor 
fluxes. It is noteworthy that since (jll)) is independent of S12, the latter does not have any 
influence on equipartition, in agreement with the findings of [20 j for the canonical case. 

If F°+F T ° = 2/3 is now assumed, then Eq.JDJ implies s 2 23 = (1-F°-2F T °)/(F°-F T °) = 
1/2, which has been first observed by [2*Tj . If, on top of this, we impose the canonical 
initial fluxes satisfying Eqs.((H), then S13 = is also uniquely required for equipartition. 
Thus, the requirement of exact equipartition of the fluxes on Earth, combined with the 
assumption that the initial fluxes are canonical, uniquely determines S23 = l/v^ and 

S13 = 0. 

The fact that the initial fluxes resulting from the astrophysical mechanism, and the 
experimentally favored neutrino mixing angles fall in the vicinity of the unique equipar- 
tition solution for the fluxes on Earth, is a striking coincidence! This is even more so, if 
we also remember the "stationary" property of these fluxes described above. 

We next turn to the fluxes as determined by the actual ranges of the mixing angles 
appearing in (jUEUHl)- Using the first order expression ©, together with (JTJ), we get 

^ = ^(1 - 2 6s 123 ) , 

F, = F T = ^(1 + 8s l23 ) , (12) 

which agrees with the conclusion of pfl |22j that for bimaximal neutrino mixing with very 
small S13, the relative neutrino flavor fluxes are F e ~ ~ F T ~ 1/3. 

Our formalism goes beyond this though, since it also considers in detail the small 
deviations from the mixing angle- values (612 — n/5.1, 623 — 7r/4, # 13 ~ 0). To linear 
order in these deviations, we find that the arriving neutrino fluxes are independent of 
<5si2, and only depend on the specific combination of <5s 2 3 and Si 3 cos<5 given in Eq.©. 
Moreover, and F T are always equal and described by the "green" line along the diagonal 
of the triangle in Fig.l. 

As seen from Egs- lfTUj) though, to second order in (<5si2, ^^23, S13), a weak intrinsic 
dependence on all four parameters (5su, ^23, S13) and cos 5 appears. To study this in 
somewhat more detail, we have reproduced in Fig.l, the implications of the constraints (@J 
|3J). The second order formulae in Eq. (jlUj) describes these constraints by the "blue-plus- 
red" region within the triangle in Fig.l, which almost completely overlaps with the "red" 

5 Barring some special cases which we have also identified, where (S12 = 0, 1/V2, 1) or (S23 = 0, 1) 
or ( s i3 — !)• Such cases are clearly excluded by experiment. In particular, the case S12 = 1/V2 studied 
in this equipartition context in has been recently reported to be excluded at the 5<r level |19j . 
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region resulted from the numerical analysis of the exact expressions |3J 122] • As seen there, 
> F T in the whole allowed region. This property can actually be shown analytically, 
both from the exact formalism or from the approximation Eq. (jlUj) which leads to 

Ffx — F T = A(1354 + 4 cos 2 5 + (5s 23 + s n cos5) 2 ) >0. (13) 

On the basis of this analysis, we find that in the canonical case 

0.28 <F e < 0.39 , 
0-30 <F M < 0.36 , 

0<AF„ T < 0.073, (14) 

and 

0.73 < *JL~(l + 3 6s lw ) < 1.21 , (15) 

which, as already said, are consistent with the results of |H I22j. 

A virtue of the present derivation, is that the effect of a future reduction of the experi- 
mental uncertainties on the mixing angles, can be straightforwardly read analytically from 
Eq.((TIlJ) or even Eq.tjHJ). As an example we note that if it turns out that e.g. 5si 23 = 0.1 
(compare Eq.(JHJ)) and that F M ~ F T , then the linear formulae (fT2J) should be adequate, 
leading to F e ~ 0.27 and F^ ~ F T ~ 0.37. In a future sufficiently large neutrino detec- 
tor such as IceCube [23], it might be possible to discriminate this case from the ideal 
prediction F e ~ = F T ~ 1/3, in the TeV-PeV energy range. 

Equal initial neutrino fluxes. Because of unitarity, if the initial relative flavors satisfy 
F e ° = F® = F®, then the final ones also obey F e = F^ = F T = 1/3, irrespective of 
the neutrino mixing angles. This is e.g. the situation roughly half a second or so after a 
supernova collapse and explosion, as suggested by typical simulations 7\. However, since 
such configurations can be simulation dependent, and keeping in mind the uncertainties 
related to supernova physics 7j, we may parameterize a deviation in the electron flavor 
from the initial flux equipartition in the form 

K = F ? = 1 -ir (16) 

The expected flux differences on Earth are then easily obtained from Eq.© as 

7 3 

F T - F e = (— - -Ss 12 - <5si 23 ) e 
16 2 

7 3 

F^-F e = (— - -5s 12 + 5si 23 ) e 

Ffj, — F T = 25s 123 e (17) 
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We call Eqs flTol IT7)1 the " Supernova-type case" , allowing it to cover also the possibility 
of TeV neutrino sources, which somehow produce roughly equal neutrino fluxes for all 
neutrino and antineutrino flavors. 

It is interesting to note that even in the linear approximation in this case, an initial \ijr 
flavor equipartition tends to be removed by neutrino oscillation effects, as can be seen from 
the third equation in Eq. (jl7|) . Moreover, a simultaneous measurement of the three flux 
differences in Eq.([17|). provided it is sufficiently accurate, would allow the reconstruction 
of the three quantities e, SS123 and Ssu- It should be clear, though, that our analysis is 
valid only if matter effects on the observed fluxes can be neglected. Obviously, this would 
not be the case for neutrinos which cross the Earth before detection [23]. Furthermore, 
star matter effects can be important in some regions of the supernova j2f 



F® = 1 case. As a last illustration we consider the rather exotic case where _F e ° = 1, 
F® = F® = 0. In the context of the linear approximation formulae Eqs.(jBJ), we get in this 

5 



case 



5s 12 



i? 3 , Ssi2 x 

F, = - + —-8s 123 , 

Fr = ^ + ^+^123 , (18) 

lb z 

where, in contrast to the previous situation, the relative neutrino fluxes have some sensi- 
tivity to 5si2 also. Using Eqs.(jH|SJ) we then find 

0.44 <F e < 0.67 , 

-0.24 <F M -F r < 0.24 , (19) 
where the uncertainties induced by Ssu and 5si23 are separated. 



To summarize, we have studied in this paper some useful analytical properties of 
the neutrino flavor fluxes, assuming just three active neutrino flavors propagating in the 
vacuum space (from the surface of the cosmic sites where they are produced, to Earth) 
oscillating among themselves with no neutrino decay processes 6 . We have first emphasized 
that for any given s\2 and cos S, the point (S23 = l/v^2 , S13 = 0) is essentially the unique 
point where (a) all three relative fluxes (F e , F^, F T ) on Earth are equal, when the initial 
fluxes are canonical, and (b) a common (si 2 - and cos 5-dependent) direction exists in the 
S23Si3-plane, [irrespectively of the values of the initial fluxes], along which all three relative 
fluxes (F e , F^, F T ) on Earth are stationary. These features have immediate consequences 
on the sensitivity to the mixing angles and flavor fluxes. Assuming then that the deviations 
of the neutrino mixing angles from their favored experimental values S12 = l/v^, S23 = 
l/y/2 and S13 = are small, we have expressed the observable neutrino fluxes on Earth 

6 The above formulae can of course be straightforwardly extended to cases including sterile neutrinos. 
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in terms of the original ones at the cosmic sites, keeping either linear or quadratic terms 
in the above angle-deviations, and assessed the validity of this approximation through 
a numerical comparison with the exact formalism. These expressions allow to extract 
information from Neutrino Astronomy data in a fairly simple way. 
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Figure 1: The relative neutrino flavor fluxes on Earth F e , F^, F T are shown in the canon- 
ical case Fg = 1/3, F® = 2/3, F® = 0. The colored regions correspond to mixing angle 
variations in the intervals 013)) using green: for the linear approximation ©; red: for the 
exact formalism result; blue plus red: for the second order expression tjlOj). The F e , F T , F^ 
coordinates of a specific point inside the triangle are obtained by reading the intersections 
with the appropriate axes, of lines emanating from this specific point along the indicated 
arrows; e.g., the illustrated point corresponds to (F e , F^, F T ) = (1/5,1/5,3/5). 
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